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Transition metal oxides (TMOs) exhibit many emergent phenomena ranging from high-temperature 
superconductivity and giant magnetoresistance to magnetism and ferroelectricity. In addition, when TMOs 
are interfaced with each other, new functionalities can arise, which are absent in individual components. 
Here, we report results from first-principles calculations on the magnetism at the BiFe03/YBa2Cu307 
interfaces. By comparing the total energy for various magnetic spin configurations inside BiFeOs, we are 
able to show that a metallic ferromagnetism is induced near the interface. We further develop an interface 
exchange-coupling model and place the extracted exchange coupling interaction strengths, from the 
first-principles calculations, into a resultant generic phase diagram. Our conclusion of interfacial 
ferromagnetism is confirmed by the presence of a hysteresis loop in field-dependent magnetization data. 
The emergence of interfacial ferromagnetism should have implications to electronic and transport 
properties. 

An interface between two different transition metal oxides (TMOs) can generate novel emergent states that 
are typically absent in its constituent TMO bulk. For example, when two TMO insulators LaAlOa and 
SrTi03 are used to form a bilayer heterostructure, a metallic state emerges at the interface. This kind of 
metallicity produced out of insulator is usually accompanied by the collective electronic phenomena, such as 
high-temperature superconductivity and colossal magnetoresistance. Therefore, the interface platform can pro- 
vide a unique opportunity for the design of interesting and controllable collective electronic properties, which 
cannot be realized with the individual constituents, and are more versatile in functionality as compared to their 
semiconductor counterparts. Since the pioneering discovery made by Ohtomo and Hwang', the research has 
immediately sparked a flurry of experimental and theoretical attempts to uncover other novel interfacial states. 
Among them, the potential to change the magnetic properties in the appropriately synthesized magnetic TMO 
heterostructures is of particular interest. Examples include the observation of a magnetic coupling between 
Lao.eySro.ssMnOs layers in a Lao.67Sro.33Mn03/YBa2Cu307_,5/Lao.67Sro.33Mn03 superlattice^ and a depression 
of the saturated magnetization in the Lao.67Cao.33Mn03/YBa2Cu307-5 superlattices' the enhancement of 
magnetoelectric (ME) coupling in some piezoelectric/ferroelectric heterostructures such as CoFe204/BaTi03'', 
and also the induction of ferromagnetism in the antiferromagnet BiFe03 (BFO) at the interface with ferromag- 
netic Lao.7Sro.3Mn03'''"'. The BF0/Lai_j:SrjcMn03 has recently been proposed as a candidate to engineer three- 
dimensional topological insulators'^. The control of interfacial ferromagnetism has also been demonstrated in 
multilayers where an insulating antiferromagnet is sandwiched between two ferromagnetic layers'". 
Technologically, this kind of research becomes especially relevant when one wants to control magnetization 
through the application of an electric field, which is mediated by ME coupling". An intrinsic ME coupling is 
expected to occur most naturally in the so-called single-phase multiferroic materials, where both the time-reversal 
and space-inversion symmetries are absent. The candidates for the multiferroic effect are BFO and RMn03 
(where R for rare earth elements). However, none of the existing single-phase multiferroic materials demonstrates 
significant and robust electric and magnetic polarizations at room temperature. In particular, materials like BFO 
and TbMn03 exhibit either commensurate or sinusoidal antiferromagnetism, this undesired property makes 
their potential technological applications limited. Therefore, the tuning of magnetism of these multiferroic 
materials when placed in contact with other TMOs is one of the most recent and fascinating research topics. 

As already mentioned above, a variety of fascinating properties have been revealed in YBa2Cu307-5-based 
junctions. It is naturally anticipated that the BFO/YBa2Cu307_<5 heterostructures could give rise to novel prop- 
erties^". In this work, we focus on the interfacial magnetism in a BFO/YBa2Cu307 (YBCO) superlattice. Within 
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Figure 1 | The crystallographic representation for the BFO/YBCO structure shown along the ac-plane (top panel) . Below it is a schematic drawing of the 
Fe interactions within BFO/YBCO superlattice structure and five representative spin configurations on the Fe atoms of BFO. The spin configuration 
AFM-GO is the same as that for the BFO bulk, for which a two-sublattice structure in the G-type AFM state is formed. For the spin configurations AFM-Gl 
through AFM-G4, the spin alignment deviates from that of the AFM-GO in the first two layers of BFO near the interface. In the first-principles simulations, 
the following lattice constants of the supercell are used: a = h = 5.564 A and c = 51.282 A. 



the first- principles density functional theory, we show that the fer- 
romagnetism is induced near the interface in BFO. We further 
develop an interface exchange-coupling model to obtain a generic 
phase diagram. By placing the magnetic exchange interaction 
strength extracted from the first-principles calculations, into the 
theoretical phase diagram, we are able to consolidate the observation 
of interfacial ferromagnetism in this composite material. 

First-principles simulations 

The ab initio calculations are performed based on the density func- 
tional theory by using the plane-wave basis set and the projector- 
augmented-wave method^' as implemented in the Vienna simulation 
package (VASP) code^^. Calculations are carried out within the local 
spin-density approximation to the exchange-correlation functional 
plus on-site Hubbard repulsion (LSDA + U) on d-orbitals of Fe. As 
in Ref 8, no repulsion is introduced for the c?-orbitals of Cu atoms 
when the YBCO is in the optimally doped regime. The band renor- 
malization effect, due to the electronic correlation in the paramag- 
netic state, should be cancelled when one looks into the relative 
energy among various spin configurations in BFO. We choose a fixed 
value of Ueff = 4 eV on Fe 3d orbitals throughout the work. A 500 eV 
energy cut-off was used to ensure the convergence of the total energy 
to 0.01 meV. For the BFO/YBCO superlattice under consideration, 
the supercell consists of 8 layers of BFO and 2 units of YBCO. In each 
unit of YBCO, there are 1 CuO chain and 2 CUO2 planes along the 
stacking direction, except that the interfacial CuO chain in YBCO is 
missing and replaced by Fe02 plane (an assumption followed from 
the scanning transmission electron microscopy analysis on LCMO/ 
YBCO superlattices^'); while the lattice constant in each plane is 
chosen such that the in-plane sublattice structure for the G-type 
antiferromagnetic (AFM) state in BFO is accommodated. The 
BrUlouin zone was sampled through a mesh of 4 X 4 X 1 fc-points. 
Due to the computational cost, all calculations are performed with 
the perfect superlattice without further atomic position relaxation of 
the structure. Since our focus in this work is on the magnetic struc- 
ture rather than the ferroelectricity, we do not expect a significant 



change from the atomic position relaxation. This assumption seems 
to be supported by the electron microscopy observation on other 
TMO interfacial materials that the atomic displacements near the 
interface are less than 0.1 A'''. Furthermore, since the BFO has a 
coUinear G-type antiferromagnetism while the normal state of 
YBCO is non- magnetic, we consider only coUinear spin polarization 
for various possible spin configurations on Fe atoms, as schematic- 
ally depicted in Fig. 1. The reason to consider various spin config- 
urations lies in the fact that the density functional theoretical 
calculations cannot automatically find out the global ground state 
when different spin states are too close in energy. The self-consist- 
ency calculations are iterated until the energy difference between two 
consecutive iterations is less than 1 X 10"^ eV. 

In Table I, we show the relative energy for the five spin configura- 
tions corresponding to those described in Fig. 1. As can be seen, the 
spin configuration AFM-Gl, for which the spins on the Fe sites are 
aligned ferromagnetically in the first layer of BFO near the interface, 
has the lowest total energy. This suggests strongly the possibility of 
ferromagnetism emerging near the BFO/YBCO interface, although 
only limited number of spin configurations are considered within the 
first-principles calculations. To prove that this spin configuration is 
indeed the globally stable state in energy, it is necessary to consider 
this spin configuration in the context of global phase diagram for the 
BFO/YBCO superlattice structure, which will be discussed immedi- 
ately below. In addition, the self-consistent results from the first- 
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Figure 2 | Magnetic moment distribution contributed from Ve 3d 
electrons in the BFO/YBCO superlattice obtained for the AFM-GO spin 
configuration (a) and the local density of states for Fe 3d orbitals in the 
first layer of BFO near the interface (solid lines) and deep into the BFO 
segment (dashed lines) obtained from the AFM-Gl spin configuration 
(b). In (a), alternating moments are indicating the G-type AFM in the two 
sublattices in each Fe02 plane. 

principles method shows only a slight reduction of Fe-3(i magnetic 
moments at the interface, as demonstrated in Fig. 2(a) for the AFM- 
GO spin configuration while the induced magnetic moment on Cu 
atoms (not shown here) in YBCO near the interface is only at the 
order of 0.01 f.iij- The slight asymmetric interfacial magnetic moment 
is due to the uniaxial shift of O atoms a little away from the Bi planes. 
These theoretical results suggest the change in magnetic properties 
occurs with one or two layers of BFO away from the interface. This 
kind of short-ranged electron evolution has been revealed by cross- 
sectional scanning tunneling microscopy (STM) on cuprate/manga- 
nite interfaces^''. Furthermore, we also show in Fig. 2(b) the local 
density of states (LDOS) for Fe 3d orbitals inside the BFO segment. 
Interestingly, although the Fe 3d states are localized deep into the 
BFO segment, those Fe 3d electrons of BFO near the interface are in 
the metallic state, which is readily accessible to an experimental test 
by the cross-sectional STM technique. By checking the electron 
charge from the first-principles simulations, we see that the total 
valence electronic charge on the Fe atoms nearest to the interface 
is larger than those that are deeper into BFO by an amount of about 
0.2. It seems to suggest the charge transfer plays an important role in 
the emergence of ferromagnetism and metallic electronic state at the 
interface. 

Effective exchange-coupling modeling 

The only slight reduction of the Fe 3d magnetic moments near the 
interface of BFO/YBCO superlattice justifies an analysis of the inter- 
facial magnetism within an effective spin-exchange model. Since the 
change in the magnetism occurs in only one or two layers of BFO 
away from the interface, the position dependent exchange couplings 
are restricted to these two layers. To investigate the various spin 
states for the BFO interface, we evaluate the spin configurations by 
evaluating the zeroth order spin exchange HamUtonian within a 
classical limit^^'^'. This gives 




Figure 3 | The 3D phase diagram for the normalized parameters based on 
the spin-exchange model for the five spin configurations of the BFO/ 

YBCO superlattice structure assuming /j' y/ | = — 1. (a) /j^^j/j' | versus 
/li/j/i I as viewed from the/2^/|/j' I =2 plane, (b) | versus JI 

as viewed from the J2 | = ~2 plane, (c) J2 | versus /i^^j/j' | as 
viewed from the /j'^j/f | =8 plane, (d) J2 y/|/i | versus /i^^j/j' | as viewed 
from the /]'y/|/i I = ~4 plane. The colors represent the borders between 

configuration regions (not the regions themselves.) The blue border 
confines the AFM-GO region from the AFM-Gl, AFM-G3, and AFM-G4 
regions. The green border confines the AFM-Gl region from all other 
regions. The yellow border confines the AFM-G2 region from the AFM- 
Gl, AFM-G3, and AFM-G4 regions. The red border confines the AFM-G4 
region from all other regions. The purple dot in panels (a) and (c) denotes 
the ground state parameters determined from the total energy calculations 
through the density functional theory. This shows that AFM-Gl is the 
ground state spin configuration. 



(1) 



where /y is the exchange parameter between spins / and j, S is the spin 
of the system, and C, — Oj is the difference between the spin orienta- 
tions. Since the spins in these systems are coUinear, cos(0; — 0^ is 
reduced to either 1 (0° or ferromagneticaUy aligned) or — 1 (180° are 
antiferromagneticaUy aligned). Table I shows the calculated classical 
energy for each spin configuration as well as the interfacial spin 
configuration. By solving the above model Hamiltonian, we plot in 
Fig. 3 to show the three-dimensional (3D) phase diagram for the spin 
configurations in the exchange parameter space as normalized by the 

parameter j/j' j . The phase diagram is represented in 3D to show the 

depth for the spin configuration regions. Although there are other 
possible spin configurations, this study focuses on the most likely 
candidates based on the total energy calculated by density functional 
theory. By fitting the classical energies for each phase with the 
density functional theory determined total energies, we can extract 
the exchange parameters /j' = -0.02302 eV, 4' =0.17329 eV, 
//- = 0.0061225 eV, and ]^ = 0.04502 eV. By placing this set of para- 
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Figure 4 | Field dependence of the magnetic moment (a) and 
magnetization (b) in the BFO/YBCO bilayer structure (red circles) as well 
as individual BFO (blue triangles) and YBCO (black squares) thin films 
deposited on a SrTiOs substrate. The data are collected at 100 K, which is 
above the superconducting transition temperature ( ~ 90 K) of YBCO 
bulk. To obtain the magnetization, as shown in panel (b) for each 
individual sample, we subtract the magnetic contribution from the SrTiOs 
substrate (see the inset to panel (a) for an example) by assuming a linear 
background throughout the entire magnetic- field range. The small 
hysteresis loop observed in the YBCO sample is due to a weak 
ferromagnetism in the SrTiOs substrate possibly produced by oxygen 
vacancies^^. 

meter values normalized to j/l'j into the phase diagram, we can 

establish the ground state phase, which is shown in Fig. 3(a) and 
(c) by the purple dot. This demonstrates that AFM-Gl is the ground 
state phase, and provides the first details that the interfacial spins in 
BFO are ferromagnetically aligned because of the electron interac- 
tions with YBCO. 

Magnetization measurement 

We have also performed the SQUID measurements on the BFO/ 
YBCO bilayer structure as well as separate BFO (—130 nm) and 
YBCO (—100 nm) thin films on a SrTi03 substrate. Figure 4 shows 
the field dependence of the magnetic moment and magnetization on 
the three structures at 100 K, at which the YBCO is in the normal 
state while the BFO bulk is already in the AFM state (Tjv - 643 K). A 
hysteresis characteristic of ferromagnetism in the BFO/YBCO 
bilayer, as shown in Fig. 4(a), is clearly seen. We note that the data 
of magnetic moment in Fig. 4(a) includes a tiny diamagnetic contri- 
bution from the SrTi03 itself, which is demonstrated representatively 
in the inset. Once we subtract this contribution within a linear back- 
ground approximation, and obtain the magnetization only for 
YBCO, BFO, and BFO/YBCO thin films, the magnetic characteristics 
as shown in Fig. 4(b) become more indicative. For our 130 nm thick 
BFO thin film, we do see a weak ferromagnetism (though the mag- 
netization data is a little noisy). This result is not inconsistent with 
the early observation that the magnetic properties in BFO thin films 



are thickness dependent^**. In contrast to the case of YBCO and BFO 
thin films, the magnetization for the BFO/YBCO bilayer structure 
exhibits a strong hysteresis loop, suggesting the ferromagnetic 
induction. Similar hysteresis loop has also been observed at 150 K 
for the BFO/YBCO bilayer structure fabricated on [(LaA103)o.3 
(Sr2TaA106)o.7] substrate^'. These experimental observations fully 
support the theoretical results obtained in the present work. 

Conclusion 

In summary, we have explored the magnetic properties of BFO/ 
YBCO superlattice structure within the density functional theoretical 
method. Our total-energy calculations together with the effective 
spin-exchange modeling have enabled to show that although the 
BFO bulk is G-type AFM, there existence of FM in the BFO near 
the interface. Our theoretical prediction has then been further sup- 
ported by the magnetization measurement in the BFO/YBCO bilayer 
structure. 
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